Force-induced desorption of a linear polymer chain adsorbed on 

an attractive surface 

P. K. Mishra, Sanjay Kumar and Yashwant Sing 
Department of Physics 
Banaras Hindu University 
Varanasi-221 005, India 

Abstract 

We consider a model of self-avoiding walk on a lattice with on-site repulsion and an attraction 
for every vertex of the walk visited on the surface to study force-induced desorption of a linear 
polymer chain adsorbed on an attractive surface and use the exact enumeration technique for 
analyzing how the critical force for desorption fc{T) depends on the temperature. The curve fc{T) 
gives the boundary separating the adsorbed phase from the desorbed phase. Our results show that 
in two dimensions where surface is a line the force fc{T) increases monotonically as temperature 
is lowered and becomes almost constant at very low temperatures. In case of three-dimensions 
we, however, find re-entrance, i. e. fc{T) goes through a maximum as temperature is lowered. 
The behaviour of the polymer chain at different values of temperature and force is examined by 
calculating the probability distribution of the height from the surface of the vertex at which external 
force is applied. 



PACS numbers: 68.18.Jk,68.47.Pe,36.20.-r 



When a long flexible polymer chain interacts with an impenetrable surface its conforma- 
tional properties are strongly modified in comparison with its bulk properties This 
is due to a subtle competition between the lowering of internal energy near an attractive 
surface and the loss of entropy due to constraints imposed by the impenetrable surface. For 
a strongly attractive surface, the polymer chain sticks to the surface, and for weak attraction 
it prefers to stay away from the surface. Thus there is a transition from the state when chain 
is mostly attached to the surface (adsorbed) to the state of detachment (desorbed) when the 
temperature is increased. The transition between theses two states is marked by a transition 
temperature Ta with the adsorbed phase for T < Ta and the desorbed phase for T > Ta. 

A model of self avoiding-walk on a lattice with on-site repulsion and attraction energy for 
every vertex of the walk visited on the surface provides an adequate model for understanding 
the adsorption-desorption transition 3, 0, 0| . We extend this model of self-avoiding walk 
to study the force-induced desorption of a linear polymer chain adsorbed on an attractive 
surface and calculate the critical force fc{T) for desorption as a function of temperature. 
Response of a polymeric chain to externally applied force can be measured experimentally 
by using techniques like optical or magnetic tweezers [6] and atomic force microscope 

We consider self-avoiding walks (SAWs) that start from a point on an impenetrable 
surface and experience a force / in a direction perpendicular to the surface at the other end. 
The applied force, because of its direction, favours desorption and one expects a critical 
force, fc{T), for desorption. At a given T when the applied force / is less than fc{T) the 
polymer will be adsorbed while for / > fc{T) the polymer will be desorbed. The curve 
fciT), therefore, gives the boundary that separates the desorbed phase from the adsorbed 
phase in the (/, T) plane. 

Let z = represents the surface and walks start from a point (origin) on the surface. In 
case of two dimensions {2-d) the surface is a line whereas in case of three dimensions {3-d) 
the surface is a plane. Since the surface is impenetrable, walks are restricted to half of the 
space {z > 0) only. We enumerated all SAWs upto a certain length on a square lattice in 
2-d and on a cubic lattice in 3-d. Let Cn{Ns, h) be the number of SAWs of vertex (a 
vertex of the walk represents a monomer of the polymer and a step of the walk, the chemical 
bond connecting the neighbouring monomers) having A^^ number of vertices on the surface, 
N — Ng vertices away from the surface and h the height of the end vertex of the walk from 
the surface. We enumerated and analyzed the series CN{Ns,h) upto A^ < 31 in 2-d and 
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< 20 in 3-d. The values of CnINs, h) found from the datas of exact enumerations for a 
given with all possible values of A^^ and h are given in table-2d.ps and table-3d.ps (see 
source file draft. tgz). 

The partition function is found from the relation 

Zn{u;, u) = Y, Cn{Ns, h)u''^u'' (1) 

where u = e'^"^''^'^ and u = e^^''^'^. Here is the energy of attraction (e^ < 0) of a vertex 
(or monomer) with the surface and / is the force acting at the last vertex of the walk in 
a direction perpendicular to the surface. In what follows, we set the Boltzmann constant 
/cb = 1 and = — 1 and express the length in unit of the lattice parameter or the length of 
a step of the walk. 

For fixed / we locate the adsorption-desorption transition temperature from the maximum 
of wS^)i=< Ns'>-< Ns >2) where < A^^ >= , N,Cr,iN,, h)u^^u^ (see Eq. 

(8)). In Fig. (1) we plot the values of g^^^^jj as a function of T{= for some values of 
/. The results given in this figure correspond to walks (a) of 31 vertices in 2-d and (6) of 20 
vertices in 3-(i. The curve ^an^'^ vs T of 2-d is found to have only one maximum for all values 
of / for which the adsorption-desorption transition occurs on varying the temperature, but 
in 3-d the curve for certain values of /(> 1) is found to have two maxima. The occurrence 
of two maxima in a ^Q^xogu) ^^^^^ indicates that the adsorption-desorption transition takes 
place at two different temperatures for the same value of /. In such a case the adsorbed 
phase is bounded by the desorbed phase from both low and high temperature sides. This is 
the case of a re-entrance. 

It is possible to obtain better estimates of phase boundaries by extrapolating for large 
A^. The reduced free energy per monomer defined as 

G{u,u)= lim ^logZAr(t^,ti) (2) 

can be estimated from the partition functions found from the data of exact enumerations 
for finite A^ by extrapolating to large A^. For N ^ oo one can, in general, write 

Zn{u,u)^ ^^{u,ufN^-^ (3) 

where ii{u,u) is the effective coordination number and 7 is the universal configurational 
exponents for walks with one end attached to the surface. The value of fj,{uj,u) can be 



estimated usiner ratio method 



with associated Neville table or any other method such as 



Pade analysis |9[ or differential approximants 2- The ratio method has recently been used 
to predict the phase diagram of a long flexible polymer chain immersed in a poor solvent 
near an attracting surface O, Isl- 

From Eqs. (2) and (3) we can write 

log fi{uj, u) = lim ^ log Zn{uj, u) = G{uj, u) (4) 

Zn{lj,u) is calculated from the data of CN{Ns,h) using Eq. (1) for given values of u and 
u ( or equivalently, T and /). From this we construct linear and quadratic extrapolants of 
the ratio of Zn{uj,u) for adjacent values of N as well as the alternate one. The values of 
fi{u,u) thus obtained for some values of force and temperature are shown in tables I and 
// for 2-d case and in table /// and IV for 3-d case. As shown in Figs (2) and (3) and from 
the tables(I-IV), the results for alternate give better convergence. We used this value 



of G{uj,u) to locate the transition temperature for given / from the maxima of Q^^^^-^i I 

d<Ns> \ 
9(log uj) ' ' 

The phase diagram thus obtained is shown in Fig. (4). We also plot in this figure the 
results found from the partition function of a chain of finite length; = 31 in 2-d and = 20 
in 3-(i. Except for small difference seen at high temperatures the agreement between these 
two values is good. We note that in case of 2-d the critical force for desorption increases 
monotonously as temperature is lowered and becomes constant at very low temperatures, 
T < 0.29. In case of 3-d we find, as mentioned above, re-entrance i. e. the critical force 
goes through a maximum as T is lowered and the adsorbed phase for / > 1 is bounded from 
both low and high temperature sides by the desorbed phase. The maximum of fc{T) occurs 
at T ~ 0.77 and its value is found to be 1.5, (see Fig. 4(6)). In both 2-d and 3-d cases the 
value of /c at T = is 1. The re-entrance in phase diagram has also been reported in 3-d in 



directed walk models 
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The nature of the phase boundary at low temperatures as shown in Fig. 4 is perhaps 
due to the specific configuration acquired by the part of chain that gets detached from the 
surface under the influence of the external force. As shown below, at low temperatures 
when the applied force is close to its critical value the detached part of the chain becomes 
a rod oriented along a particular direction. There is, therefore, no entropic contribution to 
the free energy from this part of the chain. On the other hand, at high temperatures the 
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part of the chain that gets detached under the influence of the external force remains as 
a self- avoiding walk in the half space with all lengths thermal fluctuations unless the force 
becomes sufficiently large compared to the critical force. 

To see this we calculate the probability distribution of h defined as 



for different values of / and T. The results are shown in Fig. 5 for 2-d and in Fig. 6 for 3-d 
at two temperatures. 

These figures show the difference in the behaviour of a polymer in the low and high 
temperature regions of the phase diagram. At high temperatures the thermal fluctuations 
are present in both the adsorbed and desorbed parts of the chain. As a consequence, even 
at / = when the chain is fully adsorbed, P{h) has a large width and its value at /i = is 
less than 0.5. It is only when / is sufficiently large than /c, the polymer chain gets stretched 
under the influence of force, as the force suppresses the thermal fluctuations. On the other 
hand, at low temperatures (T ~ 0.1) in both, 2-d and 3-d, the polymer chain seems to fully 
lie on the surface with /i ~ at / < /c. It is only when / is close to fc, P(h) gets broadened 
indicating that a part of the chain is pulled away from the surface. For f > fc the polymer 
chain seems to acquire shape of a rod in the direction perpendicular to the surface. At 
/ ~ /c a-U values of h appears equally probable indicating large fluctuations in the segment 
of the polymer that lies on the surface and the other segment that is away from the surface 
under the influence of the force. 

We calculated the entropy per monomer for / = using the relation 



and found (see Fig. (7)) that in 2-d the entropy goes to zero for T < 0.27 indicating that the 
chain is fully adsorbed and acquires a shape of one dimensional object. This may, however, 
not be strictly true in a long chain {N oo) due to presence of low energy long wavelength 
thermal fluctuations which are suppressed in a short chain. In 3-d entropy becomes constant 
equal to Sa — In fj, — 0.97, {fj, being the connectivity of a square lattice) for T < 0.40. This 
shows that for T < 0.4 the polymer chain is fully adsorbed on the surface and acquires a 
conflguration of a chain in two-dimensions. It therefore seems that below certain temperature 
the polymer chain is fully adsorbed and conflned to the surface. 




(5) 



s — 



d{TG{io,u^l)) 
dT 



(6) 
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From the above results one concludes that when the value of the external force / is close 
to the critical value fc a surface adsorbed polymer chain at very low temperatures (T ~ 0) 
has two parts; one is fully adsorbed on the surface or zipped with the surface and the other 
is in a shape of rod perpendicular to the surface. Thus near T ^ we can write the free 



energy as |ll|, 



12 



-Ns - TNsSa - f{N - N,) (7) 



where N — Ng represents the length h perpendicular to the surface; Sa is the entropy 
associated with the adsorbed state. Minimization of Eq. (7) with respect to A^^ gives 
fc{T) = 1 + Tsa and ^ = Sa- The value of ^ at T — is found from Fig. 4(6) to be 
~ 0.93 which is in good agreement with the value Sa = 0.97. Since in case of 2-d, Sa is zero 
as the surface is a line we get /c = 1 and ^ = in agreement with the result shown in Fig. 
4(a). 

The values of ^jf^ and calculated from relations 



<iV., > 



N Zn[uj,u)N 



^ J]iV,C^(iV„/i)cu^=ii" (8) 



and 

<h> 1 



l-Y,hCM{N„h)u''^u'' (9) 



and plotted in Fig. 8 show a sharp jump; a characteristic of a first-order transition. The 
result that the force-induced adsorption-desorption is first-order has also been found in the 
directed walk model by Orlandini et al [isj ]. 

In this paper we have used the lattice model of self-avoiding walk to calculate the values 
of critical force fc{T) for desorption of a linear polymer chain adsorbed on an attractive 
surface at different temperatures. The phase boundary that separates the adsorbed phase 
from the desorbed phase shows re-entrance in case of 3-d. It is found that for / > 1 the 
adsorbed phase is bounded from both the low and high temperature sides by the desorbed 
phase. In case of 2-d we, however, find that the critical force fc{T) attains a constant value 
for T < 0.29 and does not change on lowering the temperature. These features of the phase 
diagram can be understood from the configurations the part of the polymer chain that gets 
detached from the surface under the influence of the force acquires at the low temperatures. 
We have calculated the probability distribution of the height from the surface of the end 
vertex at which external force is applied at different values of force and temperature and 



found that for T ~ and f ^ fc the polymer chain has two parts; one is fully adsorbed on 
the surface and the other is in a shape of a rod perpendicular to the surface. 
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FIG. 1: The dependence of =< iVs^ > - < Ns on T for different forces is shown. 

Results are (a) for A?^ = 31 in 2-d and (6) for A?^ = 20 in 3-d. Note the existence of two peaks for 
/ = 1.1 in 3-d. 
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FIG. 2: The value of /XAr(a;, u) calculated from the relation /xjv(a;, u) = y ^'^ plotted 

as a function of at values of (/, r)=(0,1.3) and (0.8,0.6). The extrapolated value of pln{oj,u) 
for = has been used in Eq. (4) to calculate the reduced free energy per monomer. 
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FIG. 3: Same as in Fig. (7) but in 3-d and at values of if,T)= (0,2.45) and (1.5,0.77). 
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FIG. 4: The dependence of critical force fc{T) on T in (a) two and (6) three -dimensions. The 
star corresponds to results obtained from extrapolated values of the reduced free energy and cross 
corresponds to A?" = 20 (for 3-d) and N = 31 (in 2-d) respectively. In 3-d the re-entrance is observed 
for / > 1 and in 2-d fc{T) becomes constant equal to 1 for T < 0.29. 
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FIG. 5: The probability distribution curve for the height of the end monomer at which force is 
apphed from the surface in 2-d at T = 0.1 (a) and T = 1.0 (b). The results given in (a) are for 
N = 31 and in (6) for = 20. 
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FIG. 6: The probability distribution curve for the height of the end monomer at which force is 
applied from the surface in 3-d at T = 0.1 (a) and T = 2.0 (b). Other details are same as in Fig. 
3. 
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FIG. 7: The variation of entropy per monomer for / = in 2-d (dotted line) and in 3-c? (full line). 
The s becomes zero for T < 0.27 in 2-d and becomes nearly constant (~ 0.97) in 3-d for T < 0.4. 
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FIG. 8: The nature of variation of and has been shown with / at a fixed value of 

temperature (T = 0.1) in (a) 2-d and (6) 3-d. The results given in (a) are for N = 31 and in (6) 
= 20. 
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TABLE I: Values of |Uiv(a;, u) 
in 2-d for / = and T = 1.3 



16 2.85363 2.86096 

17 2.84831 2.83571 

18 2.84134 2.84697 

19 2.83707 2.82721 

20 2.83164 2.83607 

21 2.82815 2.82026 

22 2.82382 2.82739 

23 2.82093 2.81449 

24 2.81741 2.82034 

25 2.81499 2.80966 

26 2.81208 2.81451 

27 2.81003 2.80556 

28 2.80760 2.80965 

29 2.80585 2.80206 

30 2.80380 2.80554 

31 2.80229 2.79905 
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TABLE II: Values of iin{^^,u) 
in 2-d for / = 0.8 and T = 0.6 





IJ,n{oj,u, 


16 


5.88172 


17 


5.85734 


18 


5.83590 


19 


5.81692 


20 


5.79998 


21 


5.78478 


22 


5.77106 


23 


5.75861 


24 


5.74727 


25 


5.73690 


26 


5.72737 


27 


5.71859 


28 


5.71047 


29 


5.70294 


30 


5.69594 


31 


5.68942 



Y Zn-2{0J,u) 





5, 


.86876 


5, 


.84595 


5, 


.82587 


5, 


.80798 


5, 


.79198 


5, 


.77758 


5, 


.76454 


5, 


.75269 


5, 


.74186 


5, 


.73194 


5, 


.72281 


5, 


.71438 


5, 


.70657 


5, 


.69932 


5, 


.69257 


5, 


.68627 



Zm{u},u) 
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TABLE III: Values of HN{io,u) 
in 3-d for / = and T = 2.45 

10 5.02935 5.03992 

11 5.01239 4.98502 

12 4.99209 4.99918 

13 4.97988 4.96066 

14 4.96572 4.97078 

15 4.95643 4.94212 

16 4.94593 4.94974 

17 4.93860 4.92747 

18 4.93047 4.93346 

19 4.92453 4.91561 

20 4.91803 4.92046 

TABLE IV: Values of /Xiv(a;, u) 
in 3-d for / = 0.77 and T = 1.5 

^ u) = Miv(a;, u) = 

10 11.73428 11.662856 

11 11.59258 11.52274 

12 11.47432 11.42611 

13 11.37927 11.33262 

14 11.29805 11.26358 

15 11.23038 11.19728 

16 11.17161 11.14599 

17 11.12143 11.09692 

18 11.07728 11.05767 

19 11.03893 11.02022 

20 11.00485 10.98951 
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